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LINKING AGE, SURVIVAL AND TRANSIT TIME 
DISTRIBUTIONS 


Abstract. Although the concepts of age, survival and transit time have been 
widely used in many fields, including population dynamics, chemical engineer¬ 
ing, and hydrology, a comprehensive mathematical framework is still missing. 
Here we discuss several relationships among these quantities by starting from 
the evolution equation for the joint distribution of age and survival, from which 
the equations for age and survival time readily follow. It also becomes appar¬ 
ent how the statistical dependence between age and survival is directly related 
to either the age-dependence of the loss function or the survival-time depen¬ 
dence of the input function. The solution of the joint distribution equation also 
allows us to obtain the relationships between the age at exit (or death) and 
the survival time at input (or birth), as well as to stress the symmetries of the 
various distributions under time reversal. The transit time is then obtained 
as a sum of the age and survival time, and its properties are discussed along 
with the general relationships between their mean values. The special case of 
steady state case is analyzed in detail. Some examples, inspired by hydrologic 
applications, are presented to illustrate the theory with the specific results. 


1. Introduction 


There is considerable interest in estimating age and transit times of elements 
in a physical system or, equivalently, of individuals in a population, in disciplines 


ster 
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(Nauman 
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Preston et al.. 2000 Bongaarts and Feeney 


1925 von Foer- 


as diverse as population dynamics and demography {M’Kendrick 

2003|), ch'^ical engineering 


20081, and hydrology and geophysics (e.g., among many others 


Eriksson (: 
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(1996) 

Ginn 

(1999 
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|2006) 

Duffy 

(2010) 
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2011); 


Benettin et al. (20151). It has in fact become increasingly clear that the age, survival 


time, and the total time spent by each element in a system may provide additional 
key insights into specific aspects of a system’s behavior. This viewpoint, which 
can be considered as a time-integrated Lagrangian perspective, has been especially 


emphasized in groundwater systems ( 1 

daloszewski and Zuber 1982 

Gvetkovic and 

Dagan 

1994 Goode 1996 Ginn\ 

1999 

) and in the hydrological response of water- 

sheds 

theori 
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2006 

McDonnell and Beven. 201 
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itical I^Botter et al. 2011 van der Velde et al. 2012 Harman 
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2015 Porporato and Calabrese 2015) and field approaches \Hrachowitz et al. 

2009 

Birkel et al.\ 2011| McDonnell and Beven 2014). 


Recent discussions in the literature about the role of internal variability and ex¬ 


ternal forcing (e.g., rainfall) on the properties of age distributions (Porporato and 


Calabrese 2015), as well as the differences between age and survival time distribu¬ 


tions, their degree of statistical dependence, and their symmetry under time reversal 


(e.g., Cornaton and Perrochet (2006); Harman (2014); Benettin et al. (2015)) have 
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made evident that a comprehensive theory of age and related concepts is still miss¬ 
ing. Toward this goal, in this contribution we focus on the linkage between age and 
survival time distribution in both transient and steady-state conditions. Differently 
from what was assumed in Benettin et al. (2015), we show that age and survival 
time are in general statistically dependent quantities (the only case of independence 
being the one of time- and age-independent loss (or input) in steady state). The 
theoretical framework afforded by the evolution equation of the joint distributions 
of age and survival also provides a means to easily understand the time symme¬ 
tries between age and survival, and the derivation of the general properties of the 
transit-time distribution. 

We should warn the readers unfamiliar with the previously cited literature that, 
perhaps because of the contributions from many disciplines, the terminology which 
identifies these variables is hardly unified. For example, apart from the age, the 
definition of which seems uncontroversial (r in what follows), the survival time 
(here indicated as a) is often also indicated as life expectancy, while input and 
output rates are often also called birth and death functions. The variable with 
possibly the most appellations is the so-called transit time (T), the sum of age and 
survival time, which is also indicated as travel time, life span, total life time, and 
sojourn time. As long as the mathematical formalism is clear and the notation kept 
consistent, as we have hopefully done here, we trust that these different names will 
not confuse the readers. 

The paper is organized as follows. The evolution equation for the joint distribu¬ 
tion of age and survival time is introduced and solved in section with boundary 
conditions given by the survival time distribution at birth and the related age dis¬ 
tribution at death. The solution is used to derive the transit time distribution by a 
simple integration in section The steady state conditions are discussed in section 

Finally, we present some applications in section with the purpose of showing 
some interesting details of the theory. While most of these applications have a close 
connection to hydrological and fluid mechanic systems, they are by necessity highly 
idealized to allow us to focus on the novel theoretical results, avoiding the addi¬ 
tional complications that more realistic applications with random external forcing 
(e.g, rainfall) and spatial heterogeneities would add. 


2. Joint Distribution of Age and Survival Time 

The transit time (T) of an element of a system is the sum of the time spent 
since the entrance/birth, called the age (r), and the time that it will spend before 
exit/death, called the survival time (tr). At a given time t, each element is char¬ 
acterized by a certain age and survival (and thus transit) time, which globally can 
be described by the joint distribution (p(t, r, a). In words, ip{t, r, a)dTda represents 
the (infinitesimal) amount of elements (e.g., a mass or population number having 
age between r and r -|- dr and survival time a and tr -|- dcr at time t. 

The balance equation for the joint distribution (p(t, r, a) can be obtained con¬ 
sidering that, as the system evolves in time, (p(t, r, a) is conserved along the lines 
orthogonal to the bisector in the r, a plane, which are characterized by having 
constant T. Based on these considerations, one can readily write 

^ ^ _ cV ^ 

dt dr da 


(1) 
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Figure 1. Example of the joint distribution of age and sur¬ 
vival time in transient conditions in which the input has been 
increased in time, showing how results from a simple trans¬ 
lation of the boundary conditions. 


The equation is controlled by the boundary conditions (p{t,T,a = 0) = 
which is the survival time distribution at input/birth, and p{t,T = 0,a) = no(t, r), 
which is the age distribution at output/death. An example of the evolution of 
the joint distribution is shown in Figure [l] showing how p is simply the input 
boundary condition on the a axis, translating in time along lines of constant T 
until it crosses the r axis where cr = 0. The figure also clearly shows how the two 
boundary conditions can not be independent, as will be seen more precisely later. It 
is interesting that the contribution of input and output to the system is entirely felt 
through the boundary conditions. In more general cases, elements could also enter 
with age different from zero (immigration) or exit with a non-zero survival time 
(emigration), in which case equation 0 should also contain corresponding source 
and sink terms; these generalization however will not be pursued in this paper. 

More formally, moving along the characteristic curves, defined by ^ = 1, ^ = 
— 1 and ^ = 1, which are obviously also lines of constant transit time, it is possible 
to re-express equation Q as 


( 2 ) 


dp 

ds 


= 0 , 


so that the solution is then 


(3) p{t, r, (t) = p{t — T, T = 0, (T -I- r) = p{t + a,T + a,a = 0), 
or equivalently, 

(4) p{t, T, a) = lo{t — T, a + t) = no{t + a,T + a). 

From this it is evident how the joint distribution at a given time is simply the time 
shift of the boundary conditions and that, if time is reversed, the whole process is 
flipped, with the age playing the role of the survival time and vice versa. This type 
of time symmetry will appear frequently in the following. 
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2.1. Age and Survival Time Distributions. By integrating equation ([l|) over 


a, one obtains the M’Kendrick-von Foerster (MKVF) equation {M’Kendrick 

1925 

von Foerster 

1959 

Murray 2002 

Porporato and Calabrese 

2015), describing the 


dynamics of an age-structured population equation, 


( 5 ) 


dn{t,T) ^ dn{t,T) 
dt ^ dr 


no{t,T), 


where n{t, r) = ip{t, t, a)da is the age distribution (mass over time (age)), quan¬ 
tifying the amount of substance having age r at time t. The sink term, no(t,r), is 
the age distribution at output/death, previously introduced as a boundary condi¬ 
tion for 0 . It can be written as 


( 6 ) 


no{t,T) = fi{t,T)n{t,T), 


where fi{t, t) is the age and mass specific output rate. With initial condition n(0, t ) 
and boundary condition n{t,0) = i(t), where i(t) is the input/birth rate, the so¬ 
lution of equation 0 is \Trucco\ |1965[ \Gurtin and MacCamy \ |1974[ \Keyfitz and\ 
Keyfitz 1997) 


( 7 ) 


l{t,T) = 


l(0, T-t) e- -fo ^ ^ ^ 


_ T-) Q-!o IJ.{t-T+u,u}dv 


t > T. 


On the other hand, by integrating Q over r, a corresponding equation for the 
survival time distribution is obtained, 


( 8 ) 


dl{t,a) dl{t,a) 
dt da 


lo{t,a), 


where l{t, a) = ip{t, t , a)dT quantifies the amount of substance having survival 
a at time t. The source term Zo(f, <t) is the survival time distribution at input/birth, 
which can be expressed as 


(9) 


lo{t,a) = (3{t,a)l{t,a), 


with /3(t, a) being the survival time and mass specific birth rate. The boundary con¬ 
dition is the overall output l(t,a = 0) = o{t) and the initial condition is 1(0, a). As 
for the MKVF equation, the solution is obtained with the method of characteristics 


as 

( Z(0,t-kcr) 

t < a 

(10) 

l{t,a) = 



[ o{t + a) e- lo" 

t > a. 


By integrating again either equation 0 over T or equation Q over a, the familiar 
form of the balance equation is obtained. 


( 11 ) 

with 


dw{t) 

dt 


L{t) - o{t), 


( 12 ) 


oo oo 

'j(t) = j n(t,T)dT = j l(t,a)do 


and where the input and output can also be written as t(t) = (}{a)l{t, a)da and 

o(^) = ti-Mn(t,T)dT. 
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It should be noted that, when the solution of 0 is available, the age and survival 
distributions, n(t,T) and l(t,a), and the evolution of w{t) can be directly obtained 
by integrating the joint distribution {p{t,T,a), without need to go through the 
corresponding equations 0, 0, and ( |TI| ). Looking at these equations, it is also 
worth noting again the symmetry of the problem with respect to the time reversal, 
upon which age and survival time exchange their roles, with the output becoming 
the input and the age-specific loss function playing the part of the survival-specific 
birth function and vice versa. It is easy to see, in fact, that with these substitutions 
and t' = —t, equation ([^ and 0 are interchangeable. 


2.2. Statistical Dependence of Age and Survival. It is possible at this point 
to establish a relationship between the age-specific output and the survival-specific 
input and, in turn, discuss the conditional distributions between age and survival 
times. The latter will be essential in deriving the residence time statistics in section 
1^ To this purpose, we begin by returning to equation 0 which immediately 
furnishes the relationship between the boundary conditions by setting either r = 0 
or cr = 0 (see Figure [^, 

. . lo{t,a) = no{t + a,T = a), 

no{t,T) = lo{t-T,a = t). 


These equations, when expressed in terms of their probability density functions 
(PDFs) normalized to have area one. 


(14) 


no{t,T) , . N lo{t,cr) 

fro [t, T) = -77^ and (t, a) = 


j{t) 


L{t) 


become a relationship already obtained by Niemi (1977), 
(15) 


= oit + a)fro{t + a,T = a). 


The subscripts tq and uoi in particular, refer to the variables age at death and 
survival time at birth, respectively. They will be used explicitly when it is necessary 
to refer to them as random variables to distinguish them from age r and survival 
time cr of the entire population, as in section |4.2[ 

From ( [T^ , a relationship between the birth and loss function is then obtained 
from the definition of Uq and Iq, in equations 0 and 0, in which the respective 
solutions for n{t,T) and l{t,a) can be substituted from^^ and ([To|, giving 
(16) 

T)n(T - t) e- fo f^(u,r-t+u)du ^ g/J - I3{u,t-u)du t<T = a 


^(t, T)i(t - t) e- ff M(t-r+u,u)du ^ g- ldit-u,u)du ^ ^ ^ ^ 


This clearly shows that the age- and survival-specific birth and loss functions are 
not independent. 

Coming back to equation the joint distribution can be expressed in terms of 
the conditional distribution between age and survival time. To this purpose, it is 
more useful to consider the PDFs, instead of distributions, and define 


(17) 


fr,a{t,T,a) 


cr) 

w{t) 


fa\T{t,Cr\T) 


n{t,T) 

wit) 


fr\c,it,T\a) 


^(Co-) 

wit) ’ 


where = fritjT) and = fcrit,cr) are the marginal PDFs of age and 

survival time, respectively. 
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Focusing, as an example, on the conditional PDF of survival given age, an ex¬ 
pression can be derived combining ( [T7| ) and 

, I , _ no(t-b(T,T-f cr) 

J(t\t (^5 - 


(18) 


n(t,T) 


Substituting in (181 the solutions for no and n given by ([^ and (|^, one readily 
obtains 


(19) 


fa\r{t,a\T) = 


fj,{t + a,T + a)e- f^{u,r-t+u)du ^ ^ ^ 

fJ.{t + a,T + cr)e- t^{t-r+u,u)du ^ ^ 


Only when this expression is equal to the marginal distribution of age, are the age 
and survival time statistically independent. Thus comparing with the marginal 


( 20 ) 




/(0,1+fT) U (d(^ud+(y—u)du 
w{t) 

o(t+a) — Jg fi(t+<7—u,u)du 
w(t) 


t < a 


t > a 


one sees that, in general, marginal and conditional probability distributions are 
different. This implies that (when elements are sampled at random in the system) 
age and survival are typically statistically dependent variables. This remains true, 
in transient conditions, even when the loss and birth functions are constant, because 
the two distributions remain different for t < t and t < a. Further considerations 
on (19) will be given, for steady state conditions in section]^ 

3. Transit Time 

As already said, the transit time, T, is the total time spent by an element in the 
system, given by the sum of age and survival time, 

(21) T = T + a. 

Its distribution can thus be obtained as the distribution of the sum of the two 


random variables, the age and survival time, e.g. Springer (1979); Van Kampen 

Ht,T)= [ 

Jt 


( pOOTt , 

( 22 ) 


ip{t, r, a) dr da, 


lT<T+cr<T+dT 


which can also be written as 

T 


(23) 


(j){t,T) = j Lp{t,T',T — T')dT'= j ip(t,T — a' ,a')da'. 


0 0 
As shown in Figure]^ equation (23) is the integral along lines of constant T, which 
are orthogonal to the bisector on the plane (t, cr). Using equation ([^ one also has 
(24) 

T T T 

(j){t,T) = J ip{t + a',T,0)da'= J nQ{t + a',T)da' = J p,{t + a',T)n{t + a',T)da', 
0 0 0 
which tells us that we can know the transit time distribution at time t by summing 
up the amount leaving the system with age t = T within the time window (t, t + 
T). Only when t and cr are statistically independent is equation (24) a simple 
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Figure 2. Lines at constant transit time shown in the plane 
T,a. The boundary conditions no(t, tq) and lo{t,ao) are also 
indicated on the r and a axes, respectively. 


convolution integral. An alternative expression can be similarly obtained from 
equation ([^, 


(25) (j){t,T) = J (f{t — T + a,0,T)da= J (3{t + a\T)l{t + a ,T)da\ 

^ -T 


which instead looks back to the time window {t — r,t) and sums all the elements 
entering with survival time a = T. 

It may be useful to note that the transit time distribution refers in general to 
any element of the control volume (or population). If one instead only focuses on 
the elements entering the system (or the newborns), their transit time T is equal 
to their survival time ct, because for them r = 0, and their distribution is lo{t,a). 
Analogously, focusing on the elements leaving the system (or dying), T equals the 
age T and their residence time distribution is the equal to no(t,T). 


4. Steady State 


4.1. Distributions. Several of the previous relationships assume an interesting, 
simplified form at steady state, a necessary condition for which is that uq and Iq 
are time independent. In such a case, the balance equation for a steady state system 
is simply 


(26) 



and from equation ([^ it follows that 


(27) 


noir) = loicr 


o, 


't), 


meaning that not only the overall input equals the overall output, but also that an 
input with a fixed age r must be balanced by an output of equal survival time a. 
As a result, the joint distribution (p{T,a) is constant, and equal to no{T = T) = 
lQ{a = T) along lines of constant T. 
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In steady state conditions, age and survival time distributions are derived by 
taking t —)■ oo in 0 and ([To|, that is 


(28) 

and 


- f fj.{u)du 

n(T) = L-e ° , 


— f B(u)du 

(29) l{a) = oe o 

In particular, from ([^, the integral of (p(T,a) over t is the same as the integral 
over a, suggesting that 

(30) n(T) = l{a = t). 

In turn, it follows immediately that, at steady state, birth and loss functions must 
be equal, 

(31) ^{t) = j3{a = t). 


We note that the so-called survivor function (von Foerster 1959 Cox 1962 
Daly and Porporato. 2006|) , defined as the exceedance probability of survival at 


steady state, can be obtained by dividing either the age distribution or the survival 
distribution by the input t (or by the output o). 


(32) 


, nlr) l{a = t) - J nOdu 
S{t) = = C - L = e 0 


so that, as well known ( Cox 1962 Trucco 19651, 
dS{T 


(33) 


= froir) = = t) = fj.{T)e 0 


— f fi{u)du 


dr ^ 

The transit-time distribution at steady state can also be easily obtained by solv¬ 


ing the integral in equation (24) and substituting equation (28), 


(34) 

or, normalized as a PDF, 


— f ^{u)du 

4>{T)=Tpi{T)n{T) = iTpi{T)e » 


(35) 


/ — f a(u)du 

/T(r) =-T/i(T)e 0 

W 


Finally, regarding the conditional probabilities, from equation (18), 


(36) 
and 

(37) 


T + O- 

f . I N no(T-bcr) r , \ ~ 

/<t|t(ctt =-— =fi(T + a)e - 

n(T) 


lo(cr + T) - f PC)du 

frlairW) = -- = /3(cr -b Tjc 


l{a) 


Because of (311, the two are obviously equal. 


Comparing for example the distribution (371 with the corresponding marginal 
PDF 


( 38 ) 


L -fCu)du 

/r(r) = fa{cr = r) = -e 0 
w 
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it becomes clear that only for constant = (3, are (36) and (37) equal to their 


marginal distributions, thereby implying that age and survival time are statistically 
independent. This is essentially due to the rescaling (or memoryless) property of 
the resulting exponential distributions (Ross 2014[ ) , which is the form taken by all 
these distributions in this special case (see section 5.1). In general, however, when 
the input and loss functions depend respectively on age and survival time, the two 
variables are statistically dependent, as will be shown in detail in the applications. 


4.2. Mean values. In steady state, because of (30), the age distribution and sur¬ 
vival distribution have same mean 

(39) f = d. 

The mean age at death and mean survival time at birth are also equal, 

(40) To = CTo, 
while the mean transit time is then 

(41) f = 2f = 2a. 

With regard to the mean transit time, by definition, 


OO OO 

jT<j){T)dT jT^no{T)dT J T‘^no{T)dT 

_ O _ O _ 


(42) 


W W Tq O 

where the last equality has been obtained by multiplying and dividing by the output 
o. Now, remembering that 

OO OO 

(43) ^ J T'^no{T)dT = ^ J T'no{T)dT + var(ro) = -|- var(ro), 

O O 

where var(-) is the variance of the respective variable, then one obtains the exact 
relationship 

var(ro) 


(44) 


T = To + 


To 


Thus, in general, T > fo = (Jq, so that T represents an upper bound for both mean 
age and mean survival time. In particular, tq equals T only when the loss function 
is a Dirac delta function, for which the variance of T is zero, as in the case of a 
plug-flow system. 


In addition, substituting expression (41) into the equation above, an exact link 
between f and Tq is also obtained as 


(45) 


To , var(ro) 


2 2to 

The same condition was obtained in a somewhat different way by \Bjdrkstrdm\ 
(1978). Only in the case of /i and /3 constant, then f = tq and d = do, and 


T = 2fo = 2do, while for a Dirac delta loss function, T = tq = 

The manner in which input and loss functions depend on age and survival time 
plays a key role in determining whether the mean of the age and survival time, f 
and d, are greater or lower than the mean age at death and the mean survival time 


at birth, tq and ctq. For example, as already discussed by Bolin and Rodhe (1973) 


and Porporato and Calabrese (2015), in the case of a loss function fj, which selects 
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preferably young elements leaving older element to age in the system, the resulting 
mean age at death is lower than the mean age in the system, i.e., to<t. On the 
contrary, the case fo>T is true whenever older elements tend to be chosen by /i, 
leaving young ones to keep the mean age low in the system, compared to the mean 
age at death (see section 5.4). 


5. Applications 

We present four examples to illustrate the previously discussed theory. The first 
is a simple steady state system with constant birth and loss functions. The second 
consists of a plug-flow system in which all the elements have the same transit time. 
The third application is characterized by a periodicity of the age-independent loss 
function, while the fourth one focuses on the role of age-dependence in the loss 
function. 

5.1. Linear System at Steady State. This simplest case, which represents a 
well-mixed system, serves as a point of reference for the more complex cases pre¬ 
sented later. It is characterized by constant and equal birth and loss function. 


(46) ^ = I3 = r], 
so that the balance equation gives 

(47) w=-, 

while the normalized age and survival time distributions are 

(48) f^{T) = -qe~^'^ and f„{a) = 

with mean f = a = ^ (see Figure |3|). It is easy to show that the same exponential 


function results from (36) and 


for the conditional PDFs, so that in these 


special case age and survival are statistical independent and their joint distribution 
is simply the product of the two distributions. 

The age distribution at death (survival time distribution at birth) is simply 
obtained by multiplying n(r) by i] {l{cr) by rj), 


(49) 


noir) = T]Le 


0 1 

To = CTo = - = -• 
o 77 


which in its normalized form is equal to the age distribution, /ro(T) = rje , and 
similarly with Iq and /tq . The mean values are 

(50) 

The transit time distribution (Fig. is calculated from equation (34), 

(51) ^(T) = Tr]n{T) = 
which, written as a PDF, becomes 

(52) /t(T) = 

an Erlang-2 distribution for the sum of two independent, exponentially distributed 
random variables (see Figure]^, with mean T = ^lr\ (Cox 


1962 


Forbes et al. 


2011 ). 
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Figure 3. Age and transit time distributions for a steady 
state system with constant and age independent loss func¬ 
tion, as discussed in section The distributions are calcu¬ 
lated for r = 10 and r/ = 0 . 1 . 


5.2. Dirac Delta as Survival Time Distribution at Birth. While in the previ¬ 
ous example age and survival time were statistically independent, this application 
takes the opposite extreme of complete (i.e., deterministic) dependence. To this 
purpose, we assume that the survival distribution at birth is a Dirac delta on a 
specific survival time a*, modulated sinusoidally in time, 

(53) lo(t, (t) = A{t)d{a - a*), 

with A{t) = a + bsm{u!t). Because all the elements exit after the prescribed time 
( 7 * , this example is representative of a time-dependent plug-flow system. 

At the outlet, the age distribution at death no is the time shifted Iq, i.e., 

(54) noit,T) = A{t - t)S{t - a*), 


while inside the system, the age distribution is given by 
(55) nit, t) = (1 — 9{t — a*))A{t — t). 


where d(-) is the Heaviside function (Abramowitz and Stegun 2012|. The survival 
time distribution is 


(56) l{t, cr) = (1 — 9{<t — cr*))A{t — a* + a), 
while the joint distribution is given by 

(57) (p{t,T,a) = A{t — T)S{a + T — a*); 
see Figure]^ 

The transit time distribution can be calculated from ( [2^ , 

0 

(58) 4>{t,T) = j A(t + x)S{T — a*)dx = 

— a* 

and the conditional distribution, 

(59) cr|r) = fr\ait, T\a) = S{a + T - a*), 


5(cos(cr*w) — 1 ) 


S{T-an, 
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Figure 4. Joint distribution of age and survival time for time 
periodic plug-flow system, with boundary condition lo{t,T) 
imposed as a Dirac delta modulated by the sinusoidal ampli¬ 
tude in (53), with a = 5, 6 = 1 and a* = 10. The marginal 
distributions of age and survival time are also plotted on the 
corresponding vertical planes. 


which reflects a deterministic relationship between age and survival, imposed by 
the survival distribution at birth (see Figure]^. 


5.3. Periodic System. We now consider an extension of the first application, in 
which the input is constant but the loss function, although still independent of age, 
is now time periodic, 

(60) = r ]{ t ) = a + bsm{ujt + j ). 


From a hydrological point of view, this case is representative of a system with 
rainfall homogeneously distributed during the year but with seasonally modulated 
potential evapotranspiration and negligible other losses, 

, dw , , 

(61) — = io - n { t ) w . 


The system is still well-mixed, although the transient conditions bring about addi¬ 
tional complications that result in statistical dependence between age and survival. 

Considering t —>■ oo so that the system has forgotten the initial conditions and 
has settled on a periodic steady state, the age distribution is from Q 


(62) 


i{t,T) = toe 


- lo n(t-T+u)du _ 


b(cos(7 + u)(t-T))-coB(7 + tu;)) 


while the joint distribution can be obtained through equations 
(63) 

ip { t , r, (t) = r](t + Cr)toe“ vHdu _ ™ 


;(7 + a;(g+t))) 


The joint distribution is plotted in Figure ([^ for different days of the year 
corresponding to the four seasons. They show that the systems has low values of 
age and survival for the season with high losses (summer), while when the losses 
diminish (winter), the age and survival start increasing again. Also visible is the 
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Figure 5. Joint distributions, for a system with time peri¬ 
odic and age independent loss function ( [60| ), calculated at 
four different times of the year, t = 15 d (a), t = 105 d 
(b), t = 196 (c), t = 288 d (d). The parameters a = 0.05, 
b = 0.025, w = || and 7 = 


asymmetry with respect to the bisector indicating statistical dependence between 
age and survival induced by the time-varying conditions. The age, survival time 
and transit time distributions are plotted in Figure]^ 

The joint distribution integrated with respect to t recovers the survival time 
distribution l{t,a). However, the integral does not appear to be elementary and 
here it was only solved numerically. 

Finally, regarding the mean values (Figure [^, the loss function being age- 
independent, the mean age at death is equal to the mean age. 


(64) 


h(^)/o°°'ru(t,r)dr 

r]{t)w{t) 


The mean survival time at birth ctq was computed analytically, while the mean 
survival time a was computed numerically, and they appeared to be equal. 


5.4. Role of Age-Dependence in the Loss Function. In this last example, we 
analyze the role of the loss function under conditions of steady state. As shown 
in section]^ in steady state, the age and survival time distributions are the same, 
and therefore similar considerations also apply to a survival-time dependence of the 
birth function. We consider the age-dependent loss function 


(65) 


/r(r) = 


c -I- 1 


with c > — 1. For —l<c<0, p,isa decreasing function, thus selecting younger 
elements for output, while it is an increasing function of age for c > 0 with preference 
for older elements. 


The age distributions is obtained from (281 

(66) n(r) = 
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r, cr, T 


Figure 6. Age, survival and transit time distributions for a 
system with time periodic and age independent loss function 
( [60| ), calculated at four different times of the year, t = 15 d 
(a), t = 105 d (b), t = 196 (c), t = 288 d (d). The parameters 
are the same as in Figure]^ 



Figure 7. Time evolution of the mean age (age at death), 
mean survival time and survival time at birth and transit 
time for a system with time periodic and age independent 


loss function (60).The parameters are the same as in Figure 


whereas the age distributions at death is 
(67) no(T) = 


=-(stt)° 


For —1 < c < 0, the age PDF at death is a stretched exponential distribution 


(Sornette. 2006) and, as can be verified through (45), the mean age at death tq 


is lower than the mean age r. When c > 0, the age PDF at death is a Weibull 
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Figure 8. Age distribution, age distribution at death and 
transit time distribution for a steady state system with age 


dependent loss function given by (65) with parameter c = — \ 
(a), c = 0 (b), c = 2 (c) and c = 10 (d). 


distribution ( 

Kottegoda and Rosso 

1997 

Somette 

2006 

) and the mean age at death 

To is greater than the mean age r. In the limiting case 

c = 0, /i is a constant and 


the well-mixed system of section [SlT] is recovered. In addition, the plug-flow system 
of section [5. 2 1 is recovered when taking c —>■ oo. 

For different values of c, the age, age at death and transit time distributions are 
plotted in Figure while the joint distributions are plotted in Figure The mean 
values are shown in Figure where the mean age at death tq is lower than f when 
— 1 < c < 0, while it is larger than f when c > 0, and tends asymptotically to T 
for c —)■ 00 . As shown in equation (44), T serves as an upper bound for both mean 
age and mean age at death. 


6. Conclusions 


Our main results regard the evolution equation of the joint distribution, equation 
([^, which in turn allowed us to obtain the corresponding evolution equation for 
age (MKVF) and survival time, given by equations ([^ and (|^, respectively. The 
theory naturally led us to consider the conditional distributions of age and survival, 
equations (ffTl), (flSl), and (19), which helped us clarify some of the statements in 


the literature about the statistical dependence of these two quantities ( Cornaton 
and Perrochet\ [2006 \Benettin et al. 2015). We also obtained general relationships 
for the transit time distribution, equations (24) and (25), and discussed the sim¬ 


plifications induced by the steady state conditions; see equations (30), (27), (31), 


(35) and (37). Furthermore, we derived exact relationships among the means (44) 


and (45), although the latter was already known to Bjorkstrom (1978). 


The present theory is spatially implicit in that it considers, globally, entire popu¬ 


lations or hnite amounts of substance in a control volume ( Porporato and Calabrese 


2015). Interesting future work will consist of connecting it with the spatially explic- 


itly formulation pioneered by 

Cinn 

(1999) (and further developed by 

Weissmann 

et al. (2002) 

Cornaton and Perrochet 

(2006); Ginn et al. 

(2009 

); Cvetkovic et al. 
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Figure 9. Joint distribution for a steady state system with 
age dependent loss function /i (65), for c = 0 (a), c = 2 (b) 
and c = 10 (c). 



Figure 10. Mean transit time, mean age and mean age at 
death as function of c, given the age dependent loss function 

(p). 



amount w, that is Such dependence of birth and mortality on global 

quantities may be linked to the nonlocal nature of the pressure equation in fluid 
mechanics (Pope\ 2000 Batchelor 20001. Indeed hydrological systems are known 


to behave such that the loss function at a point is non-local (McDonnell and Beven 


2014). This is often tacitly assumed in spatially implicit, event-based formulations 


of rainfall-runoff, which use closure assumptions that depend on the total system 
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storage {Bartlett et al, 2015, in preparation). Finally, for realistic hydrologic ap¬ 
plications it is necessary to include the effect of stochasticity from the external 


environment on the input and output functions (Porporato and Calabrese 2015). 
In this case, even the mean quantities become random variables with statistical 
distributions that may be of great theoretical and practical interest. 
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